Phonon-mediated desorption of image-bound electrons from dielectric surfaces 
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A complete kinetic modeling of an ionized gas in contact with a surface requires the knowledge 
of the electron desorption time and the electron sticking coefficient. We calculate the desorption 
time for phonon-mediated desorption of an image-bound electron, as it occurs, for instance, on 
dielectric surfaces where desorption channels involving internal electronic degrees of freedom are 
closed. Because of the large depth of the polarization-induced surface potential with respect to 
the Debye energy multi-phonon processes are important. To obtain the desorption time, we use 
a quantum-kinetic rate equation for the occupancies of the bound electron surface states, taking 
two-phonon processes into account in cases where one-phonon processes yield a vanishing transition 
probability, as it is sufficient, for instance, for graphite. For an electron desorbing from a graphite 
surface at 360 K we find a desorption time of 2 • 10~^ s. We also demonstrate that depending on 
the potential depth and bound state level spacing the desorption scenario changes. In particular, we 
show that desorption via cascades over bound states dominates unless direct one-phonon transitions 
from the lowest bound state to the continuum are possible. 

PACS numbers; 52.40.Hf, 73.20.-r, 68.43.Nr 



I. INTRODUCTION 

Whenever at the surface of a solid the vacuum level 
falls inside an energy gap, that is, whenever the elec- 
tron affinity of the surface is negative, polarization- 
induced external surface states (image states) exist, as 
it is known from macroscopic electrodynamics^. Origi- 
nally predicted^, for the surfaces of liquid and solid He, 
Ne, H2, and D2 the existence of image states has by 
now been experimentally verified for a great number of 
metallio^"— as well as insulatingi^"— surfaces. In ad- 
dition, there exist a variety of dielectric materials, for 
instance, diamondi^"— , boron nitridei^, and alkali-earth 
metal oxide&iSr— , which have surfaces with a negative 
electron afhnity. They should thus support image states. 
Interesting in this respect are also electro-negative di- 
electric structures used in electron emitting devices, such 
as, cesium-doped silicon oxide Glms^^— and GaAs-based 
heterostructures2^"— . 

In contrast to intrinsic surface states^, originating ei- 
ther from the sudden disappearance of the periodic lat- 
tice potential or unsaturated bonds at the surface, image 
states are not localized at the edge but typically a few 
A in front of the solid. An external electron approach- 
ing the solid from the vacuum with a kinetic energy be- 
low the lowest unoccupied intrinsic electron state of the 
surface may thus get trapped (adsorbed) in these states 
provided it can get rid of its excess energy. Once it is 
trapped it may de-trap again (desorb) if it gains enough 
energy from the solid. Hence, in addition to elastic and 
inelastic scattering, the interaction of low-energy elec- 
trons with surfaces may encompass physisorption - the 
polarization-induced temporary binding of an electron to 
the surface. 

Unlike physisorption of neutral atoms and molecules, 
which has been studied in great detail ever since the sem- 
inal works of Lennard- Jones and collaborators^"—, ph- 



ysisorption of electrons has been hardly investigated. It 
is only until recently that we pointed out^"^ that the 
charging of surfaces in contact with an ionized g as, as 
it occurs, for instance, in the interstellar medium^"—, 
in the upper atmosphere^, in dusty laboratory plas- 
mas^i^, and in diclcctrically bounded low-temperature 
plasmas^"—, could be perhaps microscopically under- 
stood as an electronic physisorption process. 

Parameters characterizing physisorption of electrons 
at surfaces are the electron sticking coefficient So and 
the electron desorption time Tq. Little is quantitatively 
known about these parameters, although they are rather 
important for a complete kinetic description of bounded 
gas discharges (as it is in fact also the case for the sticking 
coefficient and desorption time of neutral particles which 
play a central role for the kinetic modeling of bounded 
neutral gases^"^"^^) . Very often Sc and Tc are simply used 
as adjustable parameters. 

In view of the importance of So and Te for bounded plas- 
mas, we adopted in Ref4S the quantum-kinetic approach 
originally developed for the theoretical description of ph- 
ysisorption of neutral particle a^^'^^ to calculate Sc and 
for a metallic surface. Neglecting crystal-induced surface 
states and describing the metal within the jellium model 
we obtained for an ideal surface with a classical image 
potential So ~ 10~^ and Tc ~ 10~^ s. Although So seems 



to be rather small, the product s^Tc 



10- 



which is 



the order of magnitude we expected from our study of 
charging of dust particles in low-temperature plasmas''^. 

Physisorption of an external electron implies energy 
exchange between the electron and the electronic and/or 
vibrational elementary excitations of the surface. For 
a metallic surface creation and annihilation of internal 
electron-hole pairs seem to be the main reason for elec- 
tron energy relaxation at the surface^i^. For a dielec- 
tric surface, however, the typical energy of an internal 
electron-hole pair is of the order of the energy gap, that 
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is, for the dielectrics we are interested in, a few electron 
volts. For typical surface temperatures this is way too 
large for electron-hole pairs to cause energy relaxation 
at the surface. At dielectric surfaces it has to be rather 
the creation and annihilation of phonons which leads to 
electron energy relaxation. 

For dielectrics with a large dielectric constant and a 
large energy gap, the level spacing of the two lowest states 
in the (polarization-induced) surface potential turns out 
to exceed the maximum phonon energy, which is, within 
the Debye model, the Debye energy. Hence, in contrast 
to physisorption of neutral particles, which typically in- 
volves a few bound states with energy spacings not ex- 
ceeding the Debye energji22, physisorption of electrons at 
(this type of) dielectric surfaces takes place in a deep po- 
tential supporting deep bound states whose energy spac- 
ings may be larger than the Debye energy. Relaxation 
channels involving internal electronic degrees of freedom 
being closed, because of the large gap, electron energy re- 
laxation, and hence sticking and desorption of electrons, 
has to be controlled by multi-phonon processes. 

Typical dielectric plasma boundaries are, in dusty 
plasmas^i^, graphite and melamine- formaldehyde, and 
in dielectric barrier discharges^""— Duran glass, silicon 
dioxide, and aluminum oxide. We suspect on empirical 
grounds that plasma boundaries always support image 
states, if not intrinsically then due to chemical contami- 
nation from the discharge. Based on this hypothesis we 
investigate in the following, employing a simple model 
for the polarization-induced interaction between an elec- 
tron and a dielectric surfac o^^'^^ , the desorption of an 
image-bound electron from a dielectric surface. We are 
particularly interested in how multi-phonon processes af- 
fect the competition between direct desorption, that is, 
the direct transition between bound and unbound surface 
states, and cascading desorption'^^, that is, the successive 
climbing up of the ladder of bound surface states until 
the continuum is reached. 

For the plasma boundaries just mentioned, image 
states have been so far only observed for graphitei^ (see 
Table U for the relevant material parameters). Surpris- 
ingly, the measured binding energy of the lowest im- 
age state, El^^ —0.85 eV, is lower than the energy 
of the lowest bound state in the classical image po- 
tential, which should be in fact a lower bound^. In- 
deed, for eg ~ 13.5, the dielectric constant of graphite, 
Ef « -0.63 eV. Taking either E'l''^ or Ef in conjunc- 
tion with hujo ~ 0.22 eV , the Debye energy for graphite, 
3.9 or 2.9 phonons would be required for a direct tran- 
sition to the continuum. The probability for an electron 
to de-trap from the lowest image state of graphite via 
such a transition would be accordingly small. Cascades 
using higher lying bound states as intermediaries should 
therefore become rather important. The largest energy 
difference to overcome would then be the one between the 
two lowest bound states. For the classical image poten- 
tial, e.g., this difference is 0.47 eV , implying that at most 
2.35 phonons are required for getting a cascade running 



from the lowest level. 

From these estimates we conclude that for graphite, 
with its rather high Debye energy, the number of phonons 
involved in physisorption of electrons is small enough to 
use it as an expansion parameter for the transition prob- 
ability. Taking moreover the recoil energy into account 
the dipole-active elementary excitation responsible for 
the polarization-induced surface potential imparts onto 
the external electron^, two-phonon processes even turn 
out to suffice. The other dielectrics have a much smaller 
Debye energy. The number of phonons involved is thus 
much larger. Instead of a brute force expansion other 
approaches seem to be more suitable in these cases^. 

The outline of the remaining paper is as follows. First, 
in Section II, we set up the quantum-kinetic rate equation 
for the occupancies of bound surface states and introduce 
a classification scheme for the depth of the surface poten- 
tial. In Section III we describe the microscopic model for 
the electron-surface interaction, including the static part 
which provides the surface states involved in physisorp- 
tion and the dynamic part which drives the transitions 
between these states and is thus responsible for desorp- 
tion. In Section IV we calculate the transition probability 
up to fourth order in the displacement field, thereby tak- 
ing one- and two-phonon processes into account, which 
we believe to be sufficient for graphite. This calculation 
is very lengthy''^ and cannot be totally reproduced here. 
Three appendices provide the required mathematical de- 
tails. Finally, in Section V we present and discuss our 
results before we conclude in Section VI. 



II. DESORPTION FROM MANY BOUND 
STATES 

Following Gortel, Kreuzer and Teshimai^, the kinetics 
in a manifold of bound surface states can be described by 
a quantum-kinetic rate equation. Assuming that once the 
electron is in an unbound surface state it is immediately 
pushed away from the surface, which is reasonable if we 
consider the electron as a test electron desorbing from a 
negatively charged surface (see Ref.— for more details), 
the time evolution of the occupancies of the bound sur- 
face states for the (test) electron is given by"^^ 

q' 

(1) 

where Wqq' is the probability for a transition from state 
q' to state q and Wcq = X)fe is the probability for a 
transition from the bound state q to the continuum. In 
compact matrix notation ([!]) may be rewritten as 



where n is the A^-dimensional column vector of the occu- 
pancies of the bound surface states and T is the matrix 
of the transition probabilities. 
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To determine the formal solution of this equation, 

n{t) = exp[Tt]n(0) , (3) 

the eigenvalue equation for the matrix T has to be solved. 
In general, T is not symmetric. Thus, there are right and 
left eigenvectors^, e'^ and e*^, respectively, which can be 
chosen to be orthogonal to each other. In terms of the 
right eigenvectors of T, 



(4) 



where the coefficients are determined by decomposing 
the initial distribution into eigenfunctions according to 



n(0)=5]r 



(5) 



state is at most one Debye energy beneath the continuum, 
i.e. £1 > —1, one-phonon deep if the energy difference be- 
tween the lowest two bound states is less than one Debye 
energy, i.e. A12 > —1, two-phonon deep if the energy dif- 
ference between the lowest two bound states is between 
one and two Debye energies, i.e. —1 > A12 > —2. 

Shallow and one-phonon deep potentials are typical 
for physisorption of neutral atoms and molecules. Be- 
cause of the strong polarization-induced interaction be- 
tween an external electron and a surface, physisorption of 
electrons, however, typically takes place in at least two- 
phonon deep surface potentials. Multi-phonon processes 
should thus play an important role. 



III. ELECTRON-SURFACE INTERACTION 



Due to the losses to the continuum all eigenvalues — 
turn out to be negative^. Hence, for sufficiently long 
times the image-bound electron escapes into the con- 
tinuum and the bound state occupation vanishes, i.e. 



00 



Vg. If the transitions leading to losses to 



the continuum are much slower than the transitions be- 
tween bound states, i.e. Wcq <^ Wqq' , the bound electron 
evaporates slowly into the continuum. One eigenvalue, 
Aq, is then considerably smaller than all the others and 
its right eigenvector corresponds to the equilibrium dis- 



tribution n®'^ 



,(0) 



The general solution Q can then 



be split into two terms^. 



nq{t) 



fey 



K>0 



(6) 



where the first term gives the time evolution for an equi- 
librium occupation of the bound states whereas the sec- 
ond term describes the fast equilibration of a distortion 
of the equilibrium occupation at the beginning of the 
desorption process. It is subject to much faster transi- 
tions which will be completed soon after the beginning of 
the desorption process. Since the fate of the electron for 
long times depends only on the equilibrium occupation, 
we identify the inverse of the desorption time with the 
lowest eigenvalue. 



Ao 



(7) 



This conceptual framework of desorption requires sur- 
face states and transition probabilities between them as 
input. For dielectric surfaces the transitions are driven 
by phonons whose energy scale, within the Debye model, 
is the Debye energy. It is therefore natural to measure en- 
ergies in units of the Debye energy = kBTjj. Impor- 
tant parameters characterizing the potential depth are 
then 



and Aqqf 

rtuiD 



Eg — Eql 

Hujd 



(8) 



where Eq < denotes the energy of the qth bound state. 
We call the surface potential shallow if the lowest bound 



An electron in front of a solid surface feels a 
polarization-induced attraction to the surface because of 
the coupling to dipole-active excitations of the solid. For 
a dielectric material the relevant modes are optical sur- 
face phonons^i^. 

If, for a dielectric solid with negative electron affinity, 
the kinetic energy of an external electron is less than the 
negative of the electron affinity, the electron cannot enter 
the solid, which, for the purpose of the calculation, we 
assume to fill the whole left half-space {z < 0), being ter- 
minated at 2: = with a surface whose lateral extension 
A is eventually made infinitely large. Evans and Mills- 
studied this situation by variational means. They found 
that far from the surface, the interaction potential is the 
classical image potential known from elementary electro- 
statics but close to the surface the interaction potential is 
strongly modified by the recoil energy resulting from the 
momentum transfer parallel to the surface when the elec- 
tron absorbs or emits a (dipole-active) surface phonon. 

The recoil energy makes the interaction potential not 
only nonlocal for distances less than the bulk polaron 
radius = y^h/2mujs, where m is the mass of the 
electron, eg is the static dielectric constant, and Ws = 
LUT^y (1 + es)/2 is the frequency of the surface phonon 
{ujt is the TO-phonon frequency). Most importantly, it 
makes the interaction potential finite at the surface, in 
contrast to the singular behavior of the classical image 
potential. Denoting the lateral two-dimensional momen- 
tum transfer by the simplest regular local approxima- 
tion to the true interaction potential 



Viz) 



1 TT ^ 1 

K 



-2K\z\ 



(9) 



where K is the magnitude of the vector K. It can be 
considered as a dynamically corrected classical image po- 
tential. Indeed, neglecting in the denominator the recoil 
energy, hK'^ /2m, the integral over K can be easily per- 
formed and leads to 



V,x{z) 



1 1 



A es + lz 



(10) 
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which is the classical image potential. 

The dynamically corrected image potential © is at- 
tractive. The solution of the corresponding Schrodinger 
equation will thus yield bound and unbound surface 
states. To make an analytical solution feasible, we fit 
the dynamically corrected image potential ([9]) to a 1/z 
potential that is shifted along the z axis. Forcing the two 
potentials to coincide at the surface, that is, at z = 0, we 
obtain 



V{z) 



1 1 



4 Cs + 1 z - 



(11) 



with Zc = Zg/i:. After the transformation z — > z — Zc the 
Schrodinger equation corresponding to the shifted surface 
potential reads, in dimensionless variables x = z/as and 



2An 



+ 1] (p{x) = , 



(12) 



where as = h^/me^ is the Bohr radius and Aq = 
(es — l)/4(es + 1)- Assuming that electrons cannot enter 
the dielectric surface, we solve Eq. ([T^ with the bound- 
ary condition (j){xc) = where Xc — Zc/as- The wave 
functions and energies for bound and unbound surface 
states, together with the additional boundary conditions 
we have to impose on them, are given in Appendix [Xj 

Transitions between the eigenstates are due to dynamic 
perturbations of the surface potential. The surface po- 
tential is very steep near the surface. A strong per- 
turbation arises therefore from the longitudinal acous- 
tic phonon perpendicular to the surface which causes the 
surface plane to oscillate. 

Including this type of surface vibrations and using the 
eigenstates of as a basis, the Hamiltonian for the 
surface electron can be split into three parts. 



H 



dyn 
e—ph 



(13) 



where the first term is the Hamiltonian for the electron 
in the static surface potential. 



^static ^^^^^tc^ , 



(14) 



the second term is the Hamiltonian of the free acoustic 
phonons. 



(15) 



where Q denotes a one-dimensional perpendicular wave 
vector, and the last term is the dynamic perturbation 
due to surface vibrations. Denoting for simplicity both 
bound and unbound eigenstates of the surface potential 
by \q), it is given by 



(16) 



The displacement of the surface u is related to the phonon 
creation and annihilation operators in the usual way. 



2^,ujqNs 



(17) 



with /i the mass of the unit cell of the lattice. The 
perturbation Vp{u,z) can be identified as the difference 
between the displaced shifted surface potential and the 
static shifted surface potential. Recalling PT|) and the 
transformation z —> z — Zc, it reads 



Vp{z,u) 



e^Ao , e^Ao 



z -\- u z 



(18) 



which, gearing towards a multi-phonon calculation!^, we 
expand in a Taylor series in u, 



Vp{z,u) 



-TT^--^^ +^^^ +0(u) . (19) 



IV. TRANSITION PROBABILITIES 



A. Preparatory considerations 

We intend to calculate the desorption time taking one- 
and two-phonon processes into account. Hence, we need 
to evaluate the transition probabilities Wqqi for one- and 
two-phonon processes. In general, multi-phonon pro- 
cesses have two possible origins^: (i) multi-phonon terms 
in the perturbation of the surface potential and (ii) 
multiple actions of the perturbation as it is encoded in 
the T-matrix corresponding to ^g^^/j- 

Using expansion (I19p the dynamic perturbation iJ^^^^ 



can be classified by the order in u. Up to third order, 

(20) 

where in second quantized form 

^1 = E E ^li' (bQ + b^-c?) , (21) 
Q q.q' 

V2 = - EG'.?'(Ql'Q2)(bQ,-|-bt_Qj 
<3i,Q2 q,q' 

X (bg. +b^-Q,)4c,, , (22) 

E E^?,.'(Qi'Q2,Q3)(bQ,+b^_Qj 
Qi,Q2,Q3 q,q' 

X (bQ2 + b^-Q,) (bQ3 + b^-Q,,) 4^9' • (23) 
The matrix element of the electron-phonon interaction, 

/ t \ n/2 g2jY Z^^^ 

GiAQi,...,Qn)={^] " , (24) 
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involves the electronic matrix element, 



9,9 



(25) 



whose evaluation is sketched in Appendix |Bl 

Quite generally, the transition probability from an elec- 
tronic state \q) to an electronic state \q') is given by^ 



n{q',q) 



2-K 



\{s',q'\T\s,q)t 



X S {Es - + Eg ^ Eg,) , 



(26) 



where = [ksTs)'^ with Ts the surface temperature; \s) 
and \s') are initial and final phonon states. We are only 
interested in the transition between electronic states. It 
is thus natural to average in (l26l) over all phonon states. 
The delta function guarantees energy conservation. 

The T-matrix describing the interaction between the 
external electron and the acoustic phonons obeys the op- 
erator equation. 



T — f7"iy" 



where G satisfies, 

and Go is given by 

Go = {E-H^+ie)-^ 



(27) 
(28) 
(29) 



/7" 



with Ho = atic 

For a two-phonon process we |(s', g'lTjs, q)\ in fourth 
order in u. We thus iterate T up to third order in u. 



T =Vi +V2 + V^ + ViGoVi + V1G0V2 

+ V2GaVi + ViGoViGoVi + O {u'^) , (30) 

and write for the transition probability (j26l) 
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7^(g',q) = ^7^"(<J',(J) 



(31) 



where the individual transition probabilities Tl"^{q',q) 
can be classified by their order in u. The term of O(u^), 



■.{s',q'\V^\s,q){s,q\V:W,q') , 



(32) 



gives rise to the standard golden rule approximation. 
Transition probabilities of 0{u^) vanish as the expecta- 
tion value of an odd number of phonon creation or an- 
nihilation operators is zero. We can thus drop from the 
calculation the terms 



7^l 



{s',q'\V,\s,q){s,q\V2*W,q') , (33) 

{s',q'\Vi\.S,q){s,q\V,*G*oV,*\.s\q') , (34) 

{s',q'\V2\s,q){s,q\V,*\s',q') , (35) 

{s',q'\ViGoVi\s,q){s,q\V{\s',q') . (36) 



The remaining transition probabilities are of ©(u^) and 
describe two-phonon processes. 



: 


{s\q' 


Vi\s,q)(s,q\V^\s',q') , 




(37) 


: 


{s\q' 


V^i|s,g)(s,g|K*GnT^i*|s',g') 




(38) 


: 


{s\q' 


Vi\s,q)(s,q\K*G*r,V.;\s'q') 

' i- \ 7^/ \ 7^11 U zl 




(39) 


: 


\ ; y 


yi|s,g)(s,g|yi*Gnt^i*GnFr 


s',q') , 

^ 7 y / 7 


(40) 


: 


is ,Q 


Y^ls, a) (s, alYn \s\ 0') , 

1^ 7^ / 7 




(41) 




{s',q 


'|y2|s,g)(s,ff|y*GnKls','7' 

|'.s|^7^/\ 7^l l U l 1 


) , 

/ 7 


(42) 




\* J q 










{■s\q 


'\ViGoVi\s,q){s,q\V;\.s',q') 
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(44) 




{s',q 


'\ViGoVi\s,q){s,q\V{G*oV,* 


W,q') , 


(45) 




1 1 

[s ,q 


'\ViGoV2\s,q){s,q\V,*\.s',q') 


7 


(46) 




{s',q 


'\V2GoV,\s,q){s,q\V{\s',q') 


7 


(47) 


7^" 


{s',q 


'|^iGo^iGo^i|s,g)(s,g|^i*| 




(48) 



A complete two-phonon calculation would take all 
these transition probabilities into account as they stand. 
This is however not always necessary. In the next subsec- 
tion we show that the two-phonon transition probabilities 
contain terms which are merely corrections to the one- 
phonon transition probability (j32]). Thus, for transitions 
already triggered by a one-phonon process it may in some 
cases be reasonable to neglect, in a first approximation, 
these correction terms. 



B. Calculation of the transition probabilities 

The one-phonon transition probability Ti}{q' ,q) can 
easily be brought into the form of the golden rulc'^'"' , 

2tt 



Q 

X {ub (hujQ) 5 {Eg - Eg> + Hujq) 

+ [l + nBihLjQ)]S{Eg-Eg, -hwg)} , (49) 

where the two terms in the curly brackets describe, re- 
spectively, the absorption and emission of a phonon. 

To evaluate transition probabilities numerically we as- 
sume the phonon spectrum to be adequately represented 
by the Debye model. Sums over phonon wave numbers 
can thus be transformed into integrals according to 



Q 



3Ns 



dujuj 



(50) 



Formula (jC6l) in Appendix[C]gives the one-phonon transi- 
tion probability in compact form as used in the numerical 
calculation. 

The manipulation of the two-phonon transition prob- 
abilities is rather involved and cannot be reproduced en- 
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tirely. In order to illustrate the necessary steps we take 

X (s,g|l-i*GSKk',(z') 

'K&{Es-Es- +Eq-Eq,) (51) 

as a representative example. It contains both types of 
I 



interactions: a simultaneous two-phonon interaction V2 
and two successive one-phonon interactions V\ linked by 
a virtual intermediate state arising from the iteration of 
the T-matrix. 

We begin the calculation with inserting into ([?T|) the 
expressions for Vi , V2 , and Go as given by ([2T|) , ([22l) , and 
(j29p . respectively. Inserting, furthermore, the resolution 
of the identity over electron and phonon states. 



, ^, L^s" qi,q2 Si,S2 Ql,Q2 



Es + Ea-Ho- it 



Q3 



|s2,g2) 



(52) 



Using the two identities 



(53) 



and the fact that the free resolvent is diagonal with respect to the electron-phonon states \q, s) we obtain 



n 
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('i'^1^=-T^ ^ hfj'/^ ^fjr/h e^(^-^--^)*G;^,),(Qi,Q2) 



X is e 



91 Qi,Q2,Q3,Q4 



(54) 



where all exponential factors containing electron energies have been placed in front of the phonon average. Employing 
^gjgjSst/'' — (s|e*^p'>*/'' and introducing vq — bg + bLg the above expression becomes 



ey -1 1*00 1*0 

7^'°(9','Z)=- E Y dt/he^^-^--^'^H dr/;ie^(-'-----)-G(^)^(Qi,Q2) 

X [G(}.',^(g3)G(^),(g4)]* ((t;^3(t + T)«1jJt)i;Q,^^Q.)) , (55) 



where {{■■■)) ~ X^s ^ {^\ ■ ■ - l^) / J2s" ^ average over phonon states. The operator VQ{t) evolves in 

time according to -ffp^. Hence, the four-point phonon correlation function appearing in (|55p may be rewritten as 

(56) 



and further evaluated by forming all possible contractions. Using 



(57) 
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and integrating over the times t and r finally yields 

91 Qi,Q2 



Eq - Eq^ + huQ^ - ie 
[2nB {fiujQi ) TiB {fiujQ2 ) + (^^Qi )] 



<S (£^5 - Eg, 



-2nB {fu^Qi) [1 + nB (fi^Qa)] 



(5 {Eq - + /iWQj - llhjQ^) 



[2nB {hujQ^)nB (^t^Qs) + {fnoQ.^) + (^i^^Qi 



5{Eq~Eq,) 



Eq - Eq^ - hWQ 



-2 [1 + ns {hwQj] UB (fuvQj ^ ^ 

tjq — tjq^ — hWQ^ — le 

-2 [1 + (fi^Q J] [1 + (fi^g, )] -E,-ru.^,-n.Q^, ^^^^ 

hq — tjq^ — hhjQ^ — le 



Similar expressions can be obtained for the other tran- 
sition probabilities^^. The formulae are all quite long. 
To gain more insight we classify two-phonon processes 
by the energy difference they can bridge. As can be seen 
in the above example this is controlled by delta functions 
which, quite generally, appear in two-phonon transition 
probabilities with four different arguments: S{Eq — Eqi), 
5{Eq - Eq' ± hwq), 6{Eq - Eq' ± {hujQ^ - hhjQ^)), and 
S{Eq — Eq' ± [hjjQ^ -t- fujjQ^)). For the calculation of the 
transition probabilities we can drop all contributions pro- 
portional to 5{Eq — Eqi) because in the one-dimensional 
model we are considering it implies no transition. 

In the Debye model, the maximum phonon energy is 
the Debye energy hjjjjj. Hence for terms in the transition 
probabilities that are proportional to 5{Eq — Eqi ± hujQ) 
the maximal energy difference between the initial and fi- 
nal state of the electron cannot exceed one Debye energy. 
The two-phonon transition probabilities TZ^, TZ^, TZ^, TZ7 , 
TZ}^ , TZ^^, TZ^^ and TZ^'^ have only contributions of this 
type. They are thus only corrections to the one-phonon 
transition probability TZ^ . 

Next, we consider terms proportional to S{Eq — Eq> ± 
{hujQ^ — hhjQ^)). The energies of the two phonons hojQ^ 
and fujjQ^ are both between zero and the Debye energy 
fvjj£). As they appear with different signs in the delta 
function, the energy difference between the initial and fi- 
nal state of the electron can range from —Hujd to Hlod- 
Thus these contributions do not allow to bridge levels 
that are farther apart than one Debye energy and are 
thus again merely corrections to the one-phonon transi- 
tion probability TZ^. 

Finally, we look at the contributions to the transition 
probabilities proportional to 5{Eq — Eqi±(hjjjQ^ + hjjjQr^)). 
The energy difference that can be bridged by this type 
of process is between zero and two Debye energies. Up 
to an energy difference of one Debye energy these pro- 



cesses are again corrections to the one-phonon transition 
probability. But for energy differences between one and 
two Debye energies, these are the only processes that 
contribute to the transition probability. 

This analysis leads us to divide the two-phonon pro- 
cesses into two groups, one-Debye-energy transitions and 
two-Debye-energy transitions. One-Debye-energy tran- 
sitions enable transitions between states that are at 
most one Debye energy apart, i.e. —1 < Ag^g/ < 1. 
Two-Debye-energy transitions enable transitions between 
states that are between one and two Debye energies 
apart, i.e. —2 < Ag.g' < —1 and 1 < Ag g' < 2. 

All the transition probabilities in the two phonon ap- 
proximation contribute to one-Debye-energy transitions, 
but only the transition probabilities TZ^ , TZ^'^, TZ^^ and 
TZ^"^ contribute to two-Debye-energy transitions. Among 
the contributions to the one Debye energy transitions the 
golden rule transition probability is the only "true" one- 
phonon process. 

In the following, we assume that for one-Debye-energy 
transitions the one-phonon transition probability (j49p 
dominates the corrections from the two-phonon transi- 
tion probabilities. Only when the one-phonon transi- 
tion probability is zero, which is the case for two-Debye- 
energy transitions, we will take two-phonon processes 
into account. Hence, in our numerical calculation, we use 
for one-Debye-energy transitions the one-phonon transi- 
tion probability, 

Wq'q^TZ^q'.q) , (59) 

and for two-Debye-energy transitions the two-phonon 
transition probability 

T4^g,g = TZ^q', q) + 2RcR}\q', q) + TZ'^q' , q) , (60) 

where the transition probabilities TZ? ^ TZ}^ and TZ}'^ de- 
note those parts of 7?.^, TZ^'^ and TZ^'^ which give rise to 
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two-Debye-energy transitions. They are given by equa- 
tion dCg, (lUTO)) and ([UTTj) in Appendix O The tran- 
sition probability Ti}^ does not appear exphcitly. It is 
the complex conjugate to TZ^'^ and thus subsumed in the 
second term on the rhs of Eq. ((60)) . 



be given by 



C. Regularization 

The transition probabilities for the two-phonon pro- 
cesses contain, in the present form, divergences. Specifi- 
cally within the two-Debye-energy approximation, which 
takes two-phonon processes into account only for tran- 
sitions connecting (bound and unbound) surface states 
which are between one and two Debye energies apart, 
the transition probabilities 7?.^° and T?.^"* make trouble. 

The divergences are artifacts of our one-dimensional 
model. They arise from the quantization of the elec- 
tron motion perpendicular to the surface in conjunction 
with the harmonic approximation for the lattice. The for- 
mer gives rise to arbitrarily sharp electronic energy lev- 
els for the bound states while the latter leads to infinite 
phonon lifetimes. Some divergent terms, for instance, 
/^'^2)(2, 1; 2) and I^^-^ (2, 1; 2, 2) appearing, respectively, in 

the transition probabilities 72.^°(2, 1) and 'R.^'^{2,1) (see 
Appendix [C]) , can be traced back to the diagonal matrix 
element of the linear electron-phonon interaction ((2T|) 
and could thus be eliminated with a dressing transfor- 
mation of the type used by Gortel and coworkers^. But 
other divergences, for instance, the one in the integral 
1^2) (fc'i 9; QiiQi) which appears in the rate TZ^'^{k' , q) can- 
not be removed in that manner. We decided therefore 
to regularize the divergences of the transition probabili- 
ties by taking a finite phonon lifetime into account which 
works in all cases. The drawback of this procedure is that 
it turns divergences only into resonances, whose width is 
set by the phonon lifetime, which thus becomes an im- 
portant additional material parameter. 

In order to see how a finite phonon lifetime regular- 
izes the transition probabilities, we recall bringing the 
transition probabilities into a numerically feasible form 
required to evaluate time integrals over products of time- 
dependent phonon two-point functions. For the transi- 
tion probability TZ^^ we showed this explicitly (cf. Eq. 
(IMl) and the text which followed) but the same manip- 
ulations are necessary for the other transition probabili- 
ties^. Throughout we assumed that the time evolution 
of the phonons is governed by the free phonon Hamil- 
tonian -ff^^. As a result, the two-point functions ac- 
quired an undamped time dependence. In general, how- 
ever, phonons interact, because of the anharmonicities in 
the lattice potential. A more realistic model would thus 
lead to phonon two-point functions whose time depen- 
dences are damped. Ultimately, the damping leads to 
divergence-free transition probabilities. 

To account for the damping of phonons we imagine 
the retarded and advanced phonon Green functions to 



G^^-^(Q,w) 



1 



(61) 



where jq is a decay constant arising from the phonon- 
phonon interaction and the upper (lower) sign corre- 
sponds to the retarded (advanced) Green function. Since 
the phonon four-point functions appearing in the two- 
phonon transition probabilities can be linked to these 
two functions, 7q can be incorporated into the expres- 
sions for the transition probabilities. Unfortunately, for 
the surfaces we are interested in little is known about 
the microphysics of phonons. We suggest therefore to 
use a phenomenological estimate for jq which utilizes 
material parameters which, at least in principle, could 
be measured^. 



_l_v_ 

^'^ T I nVtOD 



(62) 



with 7g the Griineisen parameter, V the volume per 
atom, and /i the shear modulus. 

In order to demonstrate how our regularization pro- 
cedure works, we consider, again as an example, a four- 
point function of the type appearing in ()55p . 



+ ((fc^ (*i)^-Q2 (i2)6Q3 (i3)&LQ, iU))) 

+ {{b~Q^ {t2)bQ, {h)blQ^ iU))) 

+ {{h-Q^{ti%,{t2)blQ^{h)bQAU))) 

+ {{b-QAti)b-Q.At2)blQ^{t3)blQ^{U))) . (63) 

First, the four-point functions have to be broken up into 
two point functions. For the first term, e.g., this means 

{{bh^itiK.it^)bQ,{ts)bQ,{U))) 

^{{bhSti)bQsit3))){{bl^Sh)bQ,{U))) 

+ {{bhAti)bQAU))){{bl^,{t2)bQAt3))) . (64) 

Because of translational invariance, the two-point func- 
tions are proportional to Sq-^^q^, even in the interacting 
case. The diagonal elements of the expectation values in 
(IMl) can be evaluated using the spectral theorem, 

{{bhit')bQm = ^ y ^ J,t,(a.)e— )d.. (65) 

and 



{{bQit)b^Qit'))) = ^ / J,t,(c.)e^"e-*-(*-*')da. (66) 



9 



with the spectral function 

«(G-«(w)-G'4(w)) 



(67) 



which contains the damping factor 7q via the retarded 
and advanced phonon Green functions. Note, we neglect 
in the spectral function the contribution proportional to 
(5(w). Since the formulae as they stand reproduce in the 
limit 7Q — >■ the correct expressions for infinite phonon 
lifetime, we conclude that the neglect of the 5{uj) term 
is in this case justified. After integration, and ([M)) 
reduce, respectively, to 



-iWQ(t-t')-7Q|t-t' 



(68) 



and 



{{bQ{t)bUt'))) = [1 + ubM] e— «(*-*')-7« 



t-t' 



(69) 



Because of the damping factors the phonon four-point 
function in (1551) becomes 



= 2nB (ujQ, )nB (u;q. )e^"«i (*+-2)-7Qi \t+r, | 



X e 



(t)-7Q2 I* 



+ 2[l + nB{iOQj] [l + nB{ujQ,)] 



X e 



-*^Ql (t+l'2)-7Qi l*+'^2 



JQ2t-7Q2 1* 



(70) 



where the ellipsis stands for terms that do not allow for 
two-Debye-energy transitions. Performing finally in (j55p 
the integral over t and t, with the phonon four-point 
function replaced by (j70p . the dominant contribution, 
that is, the term giving rise to two-Debye-energy tran- 
sitions only, can be identified as 



1 (7Qi +7Q2) 

71" (7Qi + 7Q2 )^ + K - ^q' ± i^Qi + ^Q2 ))^ 



(71) 



and, for small decay constants 7q, approximated by 



UJg, - UJq^ T 



(Wg- - ip WqJ2 + (27q^ + 7Q2)^ 



(72) 



Putting finally everything together, the corrected, 
divergence-free transition probability TZ^^ becomes 



7^^°(<7^9) = -^E E g(?!,(Qi,Qi)[gw (g2)GW,(g2) 

91 Qi,Q2 

X {2nB{hajQjnB{fiujQ2)S{Eq - Eg' + hhjQ^ + nhjQ^)g{Eq- - Eq^ ~ hwQ^,2h'yQ^ + hjQ^) 

+2[1 + nB{fiujQj][l + nB{fu^Q2)]S{Eq - Eq, -hhjQ^ -hwQ^)g{Eq, - Eq^ + nujQ^,2h^Q^ +fi-7Q2)) > (73) 

where we have used the abbreviation 



5(0,70) 



(74) 



A similar analysis can be performed for the transition probability 7?.^*. Introducing the function 



(75) 
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the corrected, divergence-free transition probability TZ is then given by 



9li92 Ql,Q2 

X {nBifiujQjnBifuxjQ2)S{Eg - Eq> + hujq^ + ftwgj 

X /(-6^92 - Eq' + SWQ2 , Eq^ - Eg' + hUQ^ , 2%Qi + fljQ^ , h-lQ^ + 2%Q2 ) 

+nB{lTi^Qi)nB{hhjQ^)5{Eq - + fi^g^ + hwQ^) 

X /(-E'92 - + , Eq, - Eg, + hwQ^ , 2h-fQ^ + h-iQ^_ , 2h-iQ^ + h-iQ^ ) 

+ [l + nB{fu^Q^)][l + nB{hujQ.J\6{Eq ~ Eg, - hwq^ - hhjQ^) 

'X-f{Eq2 - Eg' - fu^Q2,Eq^ - Eg, - nhjQ^,2h'yQ^ + h-yq^.h-fQ^ + 2h-fQ^) 
+ + nB{tu^Qi)][l + nB{hajQ2)]S{Eq - Eg, - hwq^ - hhjq^) 

^f{Eq2 - Eg, - hUJQ^,Eq^ - Eg, - fu^Q^,2hjQ^ + hjQ^,2hjQ^ + h'jQ^)) 



(76) 



TABLE I: Material parameters for graphite. 



Debye temperature 


Td 


2500if 


Dielectric constant 




13.5 


TO phonon energy 


huT 


170meV 


Gruneisen parameter 


IG 


1.7 


Shear modulus 




5 GPa 



The transition probability TZ^ has no divergence and 
hence requires no regularization. Equations (jC20p and 
(|C21[) in Appendix [C] give the final expressions for the 
corrected, divergence-free two-phonon transition proba- 
bilities as used in the numerical calculation. 



V. RESULTS 

We now insert the one- and (regularized) two-phonon 
transition probabilities into the rate equation ([1]) to 
investigate phonon- induced desorption for an electron 
image-bound to a dielectric surface. The two-phonon 
approximation is of course only applicable to dielectrics 
which have a two-phonon deep potential, for instance, 
graphite (see discussion in Section For graphite, the 
material parameters required for a calculation of the des- 
orption time are summarized in table HI Apart from the 
Debye temperature, which we varied to study the depen- 
dence of the desorption time on the potential depth, all 
numerical results were obtained for these parameters. 

A. One-phonon transition probability 

To set the stage, we start with presenting results for the 
desorption time Te calculated with one-phonon processes 
only, although these times do not apply to any of the 
dielectric materials we mentioned. 

First, we discuss the dependence of the desorption time 
on the surface temperature. If To is much larger than the 



time the electron needs to thermally equilibrate with the 
surface, prior to desorption, electron surface states are 
populated according to Ug ^ exp[— i^q/ZcTs]. Since for 
q > 2, —Eq/kTs <ti —Ei/kTs this means that the electron 
basically desorbs from the lowest surface bound state. 
Desorption requires therefore transitions from the lowest 
bound state to the upper bound states and finally to the 
continuum. They entail the absorption of a phonon. The 
likelihood of which is, according to (P^. proportional to 
the Bose-distribution ns and thus increases strongly with 
temperature. This is refiected in Fig. [T] showing that an 
increase of the surface temperature leads to an increase 
of over several orders of magnitude. Variations of the 
Debye temperature in contrast do not change the strong 
temperature dependence significantly. 

We now move on to study the effects of the poten- 
tial depth on desorption. In Section |ll] we explained 
how the potential depth can be classified by the maxi- 
mum phonon energy, the Debye energy hjj£). The relative 
depth of the potential can be changed easily by tuning 
the Debye energy whilst keeping the absolute potential 
depth constant. This is advantageous from the technical 
perspective since the cumbersome calculation of the elec- 
tronic matrix elements Z^^, does not have to be repeated. 
In order to keep the level of phonon excitation constant 
while the Debye temperature is varied we set the inverse 
surface temperature constant. 



Figure [5] shows t^^ depending on the Debye tempera- 
ture Td. The lower Td, the larger the effective potential 
depth. For a one-phonon deep potential. To > 2707 K, 
the lowest level is coupled to at least one other bound 
state by a one-phonon process. In this region the calcu- 
lation of using one-phonon processes only is applica- 
ble and leads to an increase of t^^ with increasing Td, 
that is, with decreasing effective potential depth. This is 
what one would expect as it should be easier to get the 
energy required to bridge a smaller energy difference than 
a larger one. For a shallow potential, Td > 4029 K, the 



11 



1x10 



1x10 



^ 1x10 

' CO 4 

" 1x10 



1x10' 



1x10 



1x10 



— Td 


= 3000 K 


— \ 


= 4029 K 




= 4400 K 



200 



300 



400 

T [K] 



500 



600 



FIG. 1: Inverse desorption time ^ in the one-phonon ap- 
proximation as a function of the surface temperature Ts for a 
one-phonon deep potential {To = 3000K, Td = i029K) and 
a shallow potential {Td = 44007^). 



lowest bound state can be emptied directly to the contin- 
uum. This leads to a substantial increase of t^^ signaled 
by the kink. The main conclusion at this point is that 
direct transitions which are only possible for a shallow 
potential are more effective than a cascade of transitions 
which are the only means of emptying a deep potential. 

To gain further insight, we identify the most relevant 
cascade. The one-phonon deep potential, for which the 
transition between the lowest two bound states is a one- 
phonon process, can be subdivided further, depending on 
the accessibility of the higher bound states. For the one- 
phonon deep potential the lowest level does not couple 
to the continuum directly, but it is coupled to at least 
the second bound state. Transitions between the low- 
est bound state and the third, fourth, fifth, . . . bound 
state may or may not be possible. For example, for 
Ai_2 > — 1 > Ai_3 only the second bound state, for 
Ai_3 > — 1 > Ai^4 the second and third bound state, and 
hence for Ai_„ > — 1 > Ai_„+i the second to n^^ bound 
state can be reached from the lowest state. These acces- 
sibility thresholds mark the opening of new desorption 
channels when the potential depth is reduced. Figure [3] 
shows that when the second, third, or fourth level be- 
comes available from the lowest level, t^^ increases sud- 
denly, although these steps are small. We deduce that 
the first leg of the cascade to the continuum is predomi- 
nantly the transition to the second level. 

The question arises as to how important the higher 
bound states n = 3, 4, 5 . . . are as intermediate steps for 
the second leg of a desorption cascade. To investigate 
this we calculated with different numbers of bound 
states. The image potential has an infinite sequence of 
bound states, yet for practical reasons we can take only a 
finite number of bound states into account. Figure|4]con- 
firms that adding higher bound states to the calculation 
saturates quickly. Neglecting all but a few bound 
states, say seven, is therefore justified. With two bound 
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FIG. 2: Inverse desorption time t^^ in the one-phonon ap- 
proximation as a function of the Debye temperature Td 
for different inverse temperatures 5. Tuning Td the poten- 
tial can be made shallow {Td > 4029A'), one-phonon deep 
(2707/«: <Td < 40297^) or two-phonon deep {Td < 27mK) 
for which desorption by one-phonon processes is impossible. 



states the value of r"^ amounts already to about two 
thirds of the value obtained with sixteen bound states. 
Hence the second leg of the dominant desorption channel 
is a direct transition from the second bound state into the 
continuum. The reason for the importance of the second 
level lies in the matrix element ^, which is large for 
low bound states for which more probability density is 
concentrated near the surface where the dynamic pertur- 
bation inducing desorption is strongest. 



B. Two-phonon transition probability 

Under the assumption that the true one-phonon tran- 
sition probability (j49p dominates for one-Debye-energy 
transitions the corrections coming from the two-phonon 
transition probabilities, the latter need only be consid- 
ered for two-Debye-energy transitions, for which the tran- 
sition probabilities would be zero otherwise. All the data 
presented in this subsection were obtained within this 
approximation. 

For To = 2500K, the numerical results apply to an 
electron image-bound to graphite, which has, in our nota- 
tion, a two-phonon deep surface potential. Indeed, using 
the dynamically corrected image potential, we find for 
the lowest two image states of graphite Ei = —0.347 eV 
and E2 ~ —0.114 eV. Hence, A12 = —1.06 implying 
— 1 > A12 > —2 and thus a two-phonon deep surface 
potential. Within our main assumption that the transi- 
tion probability corresponding to the minimum number 
of phonons needed to open for the first time a particular 
transition is the dominant one, the two-phonon approxi- 
mation is sufficient for graphite; n-phonon processes with 
n > 3 should yield only small corrections. 

Figure [S] concerns once more the dependence of 
on the Debye temperature, but this time also for De- 
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FIG. 3: Inverse desorption time t^^ in the one-phonon ap- 
proximation as a function of the Debye temperature for a high 
surface temperature (5 — 5). The small steps under the red 
arrows coincide with the onset of transitions from the lowest 
to the second, third, fourth, fifth, etc. bound state. 



bye temperatures leading to two-phonon deep potentials 
(Td < 2707 K). Using one-phonon transition probabili- 
ties only would drop from a finite value to zero when 
the one-phonon deep potential (A12 > —1) becomes two- 
phonon deep (A12 < —1). This happens at Tjj — 2707 K. 
Including two-phonon transition probabilities leads to a 
finite even for two-phonon deep potentials. The data 
in Fig. [5] for Td = 2500 K apply to graphite (thin ver- 
tical line). For instance, for (5 = 7, that is, Tg « 360K, 
we find « 5 • 10^ and hence a desorption time 
Te w 2 • 10-^ s. 

If a two-phonon deep potential is made shallower so 
that it becomes one-phonon deep at A12 = —1, the 
stronger one-phonon transitions set in and increases 
considerably. For high surface temperatures, for instance 
for 6 = 5, increases about fivefold. Hence, for high 
surface temperatures the one-phonon transition probabil- 
ities dominate their two-phonon corrections as expected. 
For lower surface temperatures, however, the increase 
in at the onset of one-phonon transitions becomes 
smaller. For instance, ioi S = 7 it amounts only to a 
factor of two. For very low surface temperatures, t^^ 
even drops at the threshold, e.g. for i5 = 20 by about 
40% (not shown in Fig. [SJ. In this case our assumption 
that the one-phonon transition probability dominates its 
two-phonon corrections is clearly not justified. An accu- 
rate calculation of for a one-phonon deep potential 
near A12 = — 1 requires therefore a calculation with all 
two-phonon transition probabilities included. 

Within our model for the polarization-induced surface 
potential, graphite is very close to the A12 = — 1 thresh- 
old. The neglected two-phonon corrections to the one- 
phonon transition probabilities, however, would be only 
critical if A12 were slightly larger than —1, not smaller, as 
it is in fact the case. Despite the approximations, we ex- 
pect our numerical results to be reasonable for graphite. 
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FIG. 4: Inverse desorption time t^^ in the one-phonon 
approximation as a function of the Debye temperature for 
(5 = 10 calculated with different numbers of bound states 
A^. For N = 1 desorption occurs only due to direct tran- 
sition to the continuum which dominates the rate for a shal- 
low potential {Td > 4029 K). For A'^ > 1 cascade transi- 
tions allow for desorption from a one-phonon deep potential 
(2707 K < Td < 4029 K). In this case the second bound 
state gives the most important contribution. 



especially at higher temperatures, where the resonances 
of the regularized two-phonon transition probabilities, 
which are the reminiscences of the divergences of the 
original transition probabilities, are washed out making 
the transition probabilities rather robust against small 
changes in model parameters. 

Figure [6] compares the dependence of on the sur- 
face temperature for the potential depths shallow, one- 
phonon deep, and two-phonon deep, as realized by differ- 
ent values for the Debye temperature. For all potential 
depths increases with surface temperature. The in- 
crease of TjT^ for shallow and one-phonon deep potentials 
is about the same and significantly steeper than for two- 
phonon deep potentials. For high temperatures, there- 
fore, desorption from two-phonon deep potentials lags 
behind desorption from one-phonon deep potentials. 

For the calculation with one-phonon transitions in- 
cluded only we identified the second bound state as the 
most important intermediate state for the desorption cas- 
cade. We now study the role of intermediate bound states 
when two-phonon processes are taken into account. The 
goal is again to reveal the relative importance of direct 
desorption vs. desorption via cascades. Figure [7] shows 
that saturates quickly with the number of bound 
states considered. Calculating t^^ with only the two 
lowest bound states included gives essentially the correct 
result. In the case of the one-phonon calculation, we in- 
ferred from the fact that we need only two bound states, 
that the transition from the first to the second state is the 
most important one. Within the two-phonon calculation, 
however, the interpretation is not that simple because an 
additional bound state besides the lowest one has two 
influences: Firstly, it makes cascade transitions with an 
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FIG. 5: Inverse desorption time t^^ in the two-phonon ap- 
proximation as a function of the Debye temperature To for 
a two-phonon (To < 2707 K) and a one-phonon deep poten- 
tial (Td > 2707 K) for different surface temperatures S. At 
Td = 2707 K, the onset of one-phonon transitions between 
the lowest two states, t^^ increases considerably. Data for 
Td ~ 2500 K apply to graphite (thin vertical line). 
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FIG. 6: Inverse desorption time t^^ in the two-phonon 
approximation as a function of the surface temperature 
for different potential depths. For high surface tem- 
peratures desorption from a two-phonon deep potential 
(Td = 2700 A', 2500 A', 2300 A) is significantly slower 
than desorption from a one-phonon deep potential {Td = 
2710 A, 3000 A) or a shallow potential {Td = 4100 K). 



intermediate bound state possible, secondly it alters the 
transition probability for direct transitions from the low- 
est bound state to the continuum because it also acts as 
a virtual intermediate state in the contributions to the 
transition probabilities that stem from the iteration of 
the T-matrix. Although the direct two-phonon transi- 
tion probability from the lowest bound state to the con- 
tinuum is modified it does not increase significantly with 
additional virtual intermediate bound states. Hence, the 
cascade transitions are by far more important than the 
modified direct transitions and make up almost the whole 
transition probability on their own. 

In addition to the identification of the most efiicient 
desorption channel. Fig. [7] enables us to compare our 
results with the ones obtained by Gortel, Kreuzer and 
Teshimap^. Their Figs. 1-3, 5 and 6 show for a sin- 
gle bound state as a function of the bound state energy, 
whereas our Fig. [7] shows as a function of the De- 
bye temperature, which is proportional to the inverse of 
the potential depth. Hence, apart from scaling. Fig. [7] 
is mirror inverted compared to their figures. For = 1 
our approach applies to desorption from a single bound 
state, the situation studied by Gortel and coworkers. De- 
spite the differences in the surface potential, arising from 
the fact that we are concerned with physisorption of an 
electron and Gortel et al. with physisorption of an atom, 
we also find that for potentials with depths allowing one- 
phonon transitions to the continuum desorption is much 
faster than for potentials whose depths require a two- 
phonon process for the transition from the bound state 
to the continuum. This results in a steep drop of 
at Td < 2707 K, when the one-phonon deep potential 
becomes two-phonon deep. In contrast to Gortel et al., 
however, we can include in the two-phonon calculation 



the other bound states. Then, for deep potentials, the 
stronger cascades set in and lead to a substantial increase 
ofr-i. 

Lastly we look at the relative importance of the two- 
phonon processes arising, respectively, from the expan- 
sion of the perturbation (|5D|) and the iteration of the 
T-matrix (jl9p . A two-phonon process can be simulta- 
neous, as encoded in V2, or successive, as described by 
ViGqVi. Hence, the total two-phonon transition proba- 
bility (1601) contains a contribution without virtual inter- 
mediate states, symbolically denoted by {¥2)^ (see (HH)) 
and two contributions with virtual intermediate states, 
symbolically denoted by {Vi)'^V2 and (Vi)* (see (li^ . 
dSl), and (gSl)). 

The inverse desorption time obtained from a calcula- 
tion where only two-phonon transitions due to TZ^, that 
is, due to the process {¥2)^ have been included is shown 
by the thin green and blue lines in Fig. [71 For the di- 
rect transition from the lowest bound state (A^ = 1, thin 
green line) to the continuum the process (V2)^ is always 
dominated by the processes {Vi)^V2 and (Vi)'^, as can be 
deduced by comparing the thick and thin green lines. For 
the more important cascade {N — 2, thin blue line), how- 
ever, the situation is more subtle. The processes (Vi)^ 
and {Vi)'^V2, exhibiting resonances at A12 = —1 (recall 
Eqs. ([75)1 and ([75)1 for the regularized two-phonon tran- 
sition probabilities TZ^'^ and TZ^'^, respectively) are im- 
portant only near A12 = — 1, that is, in the vicinity of 
Td = 2707 K. Far away from Td = 2707 K, it is in fact 
the process (¥2)'^ which gives the main contribution, as 
can be seen from the thin and thick blue lines in Fig. [7] 

We conclude therefore that depending on the type of 
the desorption process (direct vs. indirect via cascades), 
the energy difference between initial and final state, and 
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FIG. 7: Inverse desorption time in the two-phonon ap- 
proximation for 5 = 5 as a function of the Debye tempera- 
ture Td calculated with different numbers of bound states TV. 
Above Td ~ 2010 K the continuum is accessible from the low- 
est bound state by two-phonon processes, above Td = 4029 K 
by one-phonon processes. For Td < 2707 K the potential is 
two-phonon deep, for 2707 K < Td < 4029 K it is one- 
phonon deep and above Td = 4029 K it is shallow. For the 
thin lines labeled with (^2)^ the two-phonon process has been 
calculated using Ti? only. 



the energy of the virtual intermediate states either of 
the two-phonon processes, (1/2)^, (Vi)^V2 and (Vi)"*, may 
be the most important one and neither can thus be ne- 
glected. 



VI. CONCLUSIONS 

We investigated phonon-mediated desorption of an 
image-bound electron from dielectric surfaces using a 
quantum-kinetic rate equation for the occupancies of 
the bound surface states for the electron. To avoid 
the unphysical divergence of the classical image poten- 
tial, we included the recoil experienced by the electron 
when it couples to the dipole-active modes responsible for 
the polarization-induced interaction between the electron 
and the surface. Due to the coupling to bulk acoustic 
phonons an electron initially occupying bound surface 
states may desorb when it gains enough energy to ei- 
ther directly reach an extended state or to successively 
climb up the ladder of bound states. To allow for an 
efficient calculation of the electronic matrix elements en- 
tering the transition probabilities in the quantum-kinetic 
rate equation, we derived asymptotic approximations for 
the electron wave functions and matrix elements. 

For the dielectric materials relevant for bounded gas 
discharges (graphite, silicon oxide, aluminum oxide) or 
electron emitting devices (Cs-doped glass and GaAs het- 
erostructures) the energy spacing of at least the two low- 
est image states is larger than the Debye energy. Hence, 
provided the surface temperatures are low enough for the 




FIG. 8: Desorption channels depending on the potential 
depth. The left panel shows a shallow potential, the mid- 
dle panel a one-phonon deep potential and the right panel a 
two-phonon deep potential. For the shallow potential the pre- 
dominant desorption channel (bold orange) is a direct tran- 
sition to the continuum, for a deep potential the cascade via 
the second level is the most important channel. Full lines 
are one-phonon processes, dashed lines two-phonon processes, 
and Ed = fiLOD = ksTD- 



electron to basically desorb from the lowest bound state, 
phonon-induced desorption has to occur for these materi- 
als via multi-phonon processes, as they arise from the ex- 
pansion of the electron-surface interaction potential with 
respect to the displacement field (originating from acous- 
tic phonons) and the iteration of the T-matrix encoding 
the successive scattering of the external electron on the 
displacement field. Desorption channels involving inter- 
nal electronic degrees of freedom are closed because the 
typical surface temperatures are too low for exciting these 
high-energy modes. 

We presented results for a two-phonon deep surface po- 
tential, where the energy difference between the lowest 
two bound states is between one and two Debye energies. 
Classifying two-phonon processes by the energy differ- 
ence they allow to bridge in one- and two-Debye-energy 
transitions, we included two-phonon transition probabil- 
ities only for two-Debye-energy transitions, that is, for 
transitions where the one-phonon transition probability 
vanishes. We regularized moreover spurious singularities 
in the two-phonon transition probabilities by taking a 
finite phonon lifetime into account. 

The material parameters used for the numerical cal- 
culation apply to graphite, where the two-phonon ap- 
proximation is applicable. For a surface temperature of 
360 K, we find an electron desorption time 2-10^^ s. Be- 
sides producing an estimate for the time with which 
an image-bound electron desorbs from a graphite sur- 
face, we also investigated, as a function of the surface 
temperature and the potential depth, the relative impor- 
tance of direct vs. cascading desorption channels. For 
that purpose we used the Debye energy as an adjustable 
parameter. 

As expected, the inverse desorption time, r^T^, depends 
strongly on the surface temperature, varying several or- 
ders of magnitude when the surface temperature changes. 
Depending on the depth of the surface potential we iden- 
tified various desorption scenarios. They are summarized 
in Fig. [SI For a shallow potential, all transitions between 
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the bound states and the continuum are one-phonon pro- 
cesses. The lowest bound state can be emptied directly 
to the continuum. This is more efficient than the detour 
via higher bound states. For a one-phonon deep potential 
the transition from the lowest bound state to the contin- 
uum is a two-phonon process, whereas the transition to 
the second bound state and from there to the contin- 
uum is a one-phonon process. In this case the cascade 
with two one-phonon processes is more efficient than the 
direct two-phonon process. For a two-phonon deep po- 
tential both the direct transition from the lowest bound 
state to the continuum and the transition form the lowest 
bound state to the second bound state are two-phonon 
processes. The direct transition to the continuum is much 
slower, so that the detour via the second state is the faster 
channel. 

For most dielectrics of practical interest more than 
two phonons are required implying that for these ma- 
terials the desorption time for an image-bound electron 
may be in fact rather long. Indeed, in an ingenious ex- 
periment, using a field-effect transistor set-up with the 
gate replaced by an externally provided electron surface 
charge, Biasini and coworker a^^i^^ determined the des- 
orption time for an electron on a GaAs surface. They 
obtained 0.48 s which is rather long indeed but not un- 
expected, from our point of view, because the energy dif- 
ference between the lowest two image states of GaAs, ob- 
tained from the dynamically corrected image potential, is 
AEi2 = —0.152 eV implying more than 5 phonons to be 
necessary for that transition {huD = 0.03 eV for GaAs), 
which makes it accordingly unlikely. 

In principle, the device of Biasini and coworker a^^i^^ 
would also allow to determine the electron sticking co- 
efficient, making it a promising tool for a quantitative 
experimental investigation of physisorption of electrons 
specifically at GaAs surfaces. The empirical data about 
Te and Se are however sparse in general. In view of the 
importance these two parameters have for the complete 
kinetic modeling of bounded plasmas, it is unacceptable 
to use them as adjustable parameters, as it is currently 
the case. We encourage therefore experimental groups to 
also design devices for the investigation of electron ph- 
ysisorption at surfaces which are used or naturally appear 
as boundaries of low-temperature gas discharges. 
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Appendix A: Schrodinger equation for surface states 

To obtain bound and unbound surface states the 
Schrodinger equation has to be solved with the 

proper boundary conditions. 

First, we consider bound states for which the wave 
functions have to vanish at x — Xr and for x ^ oo. 



Substituting y = 2Aqx/k with n = ^-A^/77, Eq. p2)) 
takes the form 



K 1 

V 4 



0(y) = , 



(Al) 



whose solutions are Whittaker functions^. Hence, the 
wave functions which vanish at x = Xc and for cc — > 00 
are 



4>g{x) 



1 



where Af„ is a normalization constant defined by 



q 



dxW. 



f 2Kqx\ 



(A2) 



(A3) 



The quantum number k is determined by satisfying the 
boundary condition at the surface, W^^i/2{2AqXc/ k) = 0. 
This gives an infinite sequence of real numbers n whose 
differences are however roughly one so that we can map 
them one-to-one onto integers q, that is, q -H- k. The 
energy of the bound state labelled with q is thus given 
by 



me'^Ap 1 
2ft2 



(A4) 



Continuum states vanish only at a; = a:c. Since for 
them the energy is positive k is imaginary. It is thus 
more convenient to label them with the real number k — 
The wave function for continuum states is given 
by a linear combination of Whittaker functions. 



1 



-_{2iAokx) +cW± i(-2iAofca;) 



(A5) 



where the constant. 



c = -W_:^ l{2iAokxc)/W^_ l{-2iAokxc) , (A6) 

is chosen to enforce the boundary condition (f>{xc) = 0. 
Normalizing the wave function in a box of length L leads 
to the normalization constant A/fe = \/2Le^. The en- 
ergy of the continuum states, finally, is given by 



Eh — 



^0 1,2 



2ft2 



(A7) 



In the limit L 00 sums over continuum states can thus 
be transformed into integrals according to^ 



fc>0 



AqL 



dfc. 



(A8) 



where the L factor in front of the integral cancels with 
the ^/L factor contained in the normalization constant 
. 
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Appendix B: Evaluation of the matrix elements 

The electronic matrix element qi for two bound 
states labelled, respectively, by q and g', reads 

1 r . j^2Ao5' 



dxW^ 



X —W^, 



(Bl) 



For the efficient numerical evaluation of Afq and Z^j-^, we 
use an expansion of Whittaker functions for k > in 
terms of Laguerre polynomials^^, 



k{k — l)e 2^Ln{x) 
n'^o ii^-n)iK-n- l)r(2 - k) 



(B2) 



To calculate the bound state matrix element Z1^^, we 
can compute one matrix element after another. First, 
we choose two states q and q' which fixes the quantum 
numbers k and k' in the bound state wave functions. 
Then, we integrate over x and obtain plain numbers for 
the matrix element which can be directly inserted into 
the calculation of the transition probabilities. 

The evaluation of the electronic matrix element is more 
demanding if either one or both states are continuum 
states. A wave function in the continuum is labeled by a 
real number k. If we were only interested in the value of 
the matrix element for some k we could follow the same 
strategy as for bound states. Some transition probabili- 
ties however contain sums over all electronic states, e.g., 
the sum over qi in TZ^'^ (see Eq. (|58p ). which for con- 
tinuum states implies an integral over k. Hence k is not 
merely a parameter that we can specify in advance. It is 
rather a variable in a matrix element which thus becomes 
a function of fc, 



^q,k ~^ Zq{k) 



1 

dx<f>q{x) ^ Mx) ■ (B3) 

X" 



For the two-phonon calculation we specifically need the 
bound-continuum matrix element Z^'^, for a given q and 
< k < oo and the continuum-continuum matrix ele- 
ment Zl f^, for < fc < oo and k' small. Because of the 
complicated structure of the Whittaker function the ma- 
trix elements cannot be obtained straightforwardly. To 
make their calculation feasible we constructed approxi- 
mate expressions for the continuum wave function (j)^. (x) 
for small and large fc, respectively, calculated the ma- 
trix elements in these two limits, and then interpolated 
between them with a Fade approximation^^. 

We begin with the limit fc — > 0. In this limit, the 
Schrodinger equation for the continuum states is 



2 A 

(l)"{x) + —(t>{x)^0 , 



(B4) 



which after the substitutions t = 2v^2Ao^ £^nd (p = tip 
takes the form of the Bessel differential equation^. 



r 



(B5) 



Hence, in the limit fc — > 0, the continuum wave func- 
tion satisfying the boundary condition (j)k{xc) = can 
be written as a linear combination of Bessel and Neu- 
mann functions. 

There is however one technical caveat. Bessel and Neu- 
mann functions grow beyond limit for large x (see formula 
(9.2.1) in ReiM) and cannot be normalized. Within the 
matrix element this is not dangerous because the decreas- 
ing factors 1/a;" compensate the divergence at large x. 
The only problem left is to find the prefactor by which 
we have to multiply the linear combination so that it has 
for small x the amplitude of the correct (jjk^o (x) ■ Once 
we have this factor, the linear combination is normalized 
in the sense that its envelope coincides for small x with 
the envelope of the correct (/)fe_i.o(x). 

The most direct way to obtain the required multiplica- 
tion factor is to perform the limit fc — )• in the continuum 
wave function (jA5|) . Due to the complicated structure of 
the Whittaker function this is however not feasible. In- 
stead, it is better to determine the factor from the solu- 
tion of Eq. ((B5|) . 



(x)fc^o = 2v/2AoxJi(2v/2Aox) 



(B6) 



which satisfies the boundary condition 0(0) = 0. The 
solution of ([T2|) satisfying the same boundary condition 
is given by 



(t)k{x) = A/jT^Af.^ i(2iAofcx) 
with the normalization constant 



^/ Lfc(l-e-2-A) 



(B7) 



(B8) 



In the limit fc (|B7p merges into (jB6|) because 



M 



(2iAofcx) !■ iky/2AoxJi{y/8AQx) . (B9) 



fe^O 



For the normalization constant we thus obtain A/'j. ^ — > 
^/tt/ [Lk] for fc 0. The fact that fc and x are tied 
together in the argument of M_j/j, i/2(2iAofcx) suggests 
that an approximation for small fc is at the same time 
an approximation for small x. Since we need the wave 
function only for small x and small fc, we expect the 
replacement (|B9P to provide good results for the desired 
normalization factor. 

Since the amplitudes of Bessel and Neumann functions 
are the same for large x, we can use the normalization 
obtained for the Bessel function for the Neumann func- 
tion as well. A normalized approximation satisfying the 
correct boundary condition, (pkixc) — 0, is then given by 



X Ji(2v/2Aox) - 5iVi(2v/2Aoa;) (BIO) 



with 



-Ji(2V2AoXe)/iVi(2V2AoXc) 



(Bll) 
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Using (jBlOl) the matrix element Z^j. can be evaluated 
in the limit A; — >■ 0. We find 



dx 



(j)q{x)(l)k^o{x) 



ri-l/2 



, (B12) 



where the k dependency is separated so that the remain- 
ing integral gives a fc-independent quantity, a^, which 
has to be obtained numerically. Similarly, 



= dx 



4>k' ^a{x)4>k^a{x) 



'kk' 



aBX^ 



LCLB 



(B13) 



with again a constant to be determined numerically. 

We now proceed to the approximation of the contin- 
uum wave functions for large k. The higher the energy 
of the continuum states, the lesser the potential at the 
surface changes the plane wave behavior far from the sur- 
face. Therefore, in the limit fc — oo we can accurately 
describe the wave function by a plane wave. Using the 
asymptotic form for the Whittaker function (page 116 in 
Ref^), 



(B14) 



the continuum wave function in the limit k ^ oo satisfy- 
ing the correct boundary condition, 4>k{xc) — 0, can be 
approximated by 



'2/L sin(Aofc(a; — Xc)) 



(B15) 



Employing, finally, the Fourier integral^ the matrix 
elements for large k can be shown to be^ 



^q.k^ca 



and 
72 



ji-1/2 



\fLMq (Aofc)3a; 

6-4 



-ttM^' 1 

n+1 K,-^ 



2Ao 



(BI6) 



1 



32 



Las V 1 + 5' k^'J^xl 



(B17) 



independent coefficients and 



have to 



where the k - 

be worked out again numerically. 

Having calculated the leading terms for the matrix el- 
ements for small and large k, we can combine these two 
limits via a Fade approximation in terms of ^/k. The co- 
efficients are chosen in such a way that the Fade approx- 
imation matches the leading term of both limits: fc 
and k ^ 00. Then, the matrix elements read 



1 



yn 



where /3" = ct'^/b'i, and 



7^ — 



1 

LaB 1 + /32F/2 ' 



(B18) 



(B19) 



where Pi = a1/b1 



Appendix C: Transition probabilities in a compact 
form 



In this appendix we list the one- and two-phonon tran- 
sition probabilities as used in the numerical calculation 
of the dcsorption time. We implicitly assume that q la- 
bels both bound and unbound states. For a bound state 
q is simply an integer (to be mapped onto «(?)) whereas 
for a continuum state q stands for a real number k. To 
obtain a compact form for the transition probabilities it 
is moreover convenient to introduce dimensionless vari- 
ables, 



E„ 



D 



v{x) ■ 



hijj 

Eq — Eqi 

k^ ' 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 



The transition probability of 0{u), that is, the one- 
phonon transition probability employed in a golden rule 
approximation, is given by 



(9, 9) = ^ E (Q) (Q)] * ("s (^q) ^ (^9 - ^1' + ^q) + [1 + (^q)] ^ iEq - Eq, - hWQ)) 



3^2^72 72 



(C6) 
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where 



e-*^-^' - 1 



for < -Ag^g, < 1 , 



and 



9.9' 



for < Ag^g, < 1 . 



(C7) 



(C8) 



Otherwise I^^^ and l'^^^ are zero. Depending on whether q and q' denote bound or continuum states, the electronic 
matrix element ^, is either given by (IB1|) . (jBlSp . or (|B19|) . 

As explained in Subsection IIVBI we keep only those parts of the transition probabilities of 0{u*) which give rise to 
what we call two-Debye-energy transitions, which are transitions between states that are between one and two Debye 
energies apart. In the unregularized form, that is, in the form which diverges in particular situations (see Subsection 
IIV CI for a discussion), the two-phonon transition probabilities included in our calculation are given by 

{q': 9) = IT E ^9',9 (Qi^Qi) (^2, Q2)] * (2^5 {fujQjnB {f^Q,) S {Eg - Eg, + hujQ, + hWQ,) 



-2 [1 + UB (^WQi)] [1 + ns (^WQs)] ^ [Eq - Eg, - hojQ^ - hujQ^)) 



7^l" iq',q) 



z. 



(C9) 



G\,^g{Q^Mx){G\^g{Q-,)G\,^g^{Q^)X 

S {Eg - Eg, + hujQ^ + hojQ^ 



'2 [1 + hb (fiujQ^ )] [1 + hb (hujQ 



Eg - Eq^ + hwQ^ - ie 

S {Eg - Eg, - hhjQ^ - huQ^) 



Eg Eg-^ hhJ, 



^ D r/1 



(CIO) 



n'\q',q)^^Yl E G\,^g,{Q^)Gl^g{Q^)[G\,^g^{Q2)Gl^g{Q2)\' 



91,92 Qi,Q2 
X ( ns (fiwQ J ns (?icjQ J 



5 {Eg - Eg, + hwQ^ + hiOQ^) 



+nB {hajQ^)nB {fiojQ^) 



[1 + riB (ftwgj] [1 + UB (ftwQa)] 
[1 + TIB (ftwgj] [1 + UB (ftwQa)] 



{Eg - Eg^ + hhjQ^ + ie){Eg - Eg^ + Huj Q ^ - ie) 

S {Eg - Eg, + hUQ^ + hUjQ^) 

{Eg - Eq^ + SwQi + i€){Eg ~ Eg^ + ftwQj - ie) 

5 {Eg - Eg, - hOjQ^ - hUQ^) 



{Eq - Eq^ - hhjQ^ + ie){Eg ~ Eg^ - hwq^ - ie) 

5 {Eg - Eg, - hwQ^ - hojQ^) 

{Eq - Eq, ~ hujQ^ + it){Eg - Eg^ - hojQ^ - ie) 



E 4,91 4,94,92 42.9' {^{2) 9; 91, 92) + /(% 9; 91, 92) + I^) (9', 9; 91, 92) 



+-ff2) (9',9;9i,92)) , 
where, in the limit e — >■ 0, the auxiliary integrals are defined by 



^(2) (9', 9) = 



dx- 



.5x 



Ja ,-i 



le*(-N,,'-) -1 
1 



oSx 



1 



(A, 



9.9 



X 



1 + 



1 

=^(A,,,'--) 



(Cll) 

(C12) 
(C13) 
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1(2) W^Q'^Qi) 



1(2) ii'^Q'^qi) 



dx- 



dxx 



1 eK~'^i.i'-^) - 1 \,qi +x-ie 



o6x 2 



A, 



9,91 



(C14) 
(C15) 



I{2) (9',?; 91, 92) 



I{2) (9','?; 91, 92) 



dx 



g<5x _ I ^5(-A, _ ^ Ag>,q^ - X + ie Ag^q^ + x - ie 



dx- 



A ,_i e*^ - 1 g5(-A, _ ^ Aq^q^ +x + ie Aq^q^ + X - le 



(C16) 
(C17) 



^(2) (9, 9; 91, 92) 



^f2) (9', 9; 91, 92) 



dxa; 



1 + 



1 + 



1 



1 



Aq',gi + X + ie Ag^g2 — a; — «e 



a; + ie Aq^q^ — x — ie 



(C18) 



(C19) 



The two-phonon transition probabilities TZ^^ and TZ^'^ diverge for the situations studied in this paper. They have to 
be therefore regularized. The regularization procedure outlined in Subsection IIV Cl leads then to the divergence- free, 
regularized transition probabilities, 



7^^"(9',9) 



4.94.91^91.9 (A'2)(9''9;9i) +/('2)(9',9;9i)) 



(C20) 



^''^'^'''^^ ^ 2hW E 4.91 4.94.9. 4.9' (^(2) (9',9;9i,92) + (9', 9; 9l , 92)) , (C21) 

D rn 



where the integrals are given by 
/■I 

r3' t' 



I(2)i<l ,9; 91) 



1(2) (9', 9; 91) 



dx- 



x 



A„ „, -1 



dxx 



.5x 



91.92 



^ g{Aq^q,+X,v{-X- Aq^q^)+2U{X)) , 



1 



9{^q,qi - a;, i^(Ag^g/ - x) + 2v{x)) , 



(C22) 
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1(2) (9', 9; 91, 92) 
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5x 



le*(-A,,,'--) _1 

X [/(Ag',qi - X, Aq^q^ + X, Vix) + 2v{-X - Aq^q,), 2v{x) + v{-X - Aq,q,)) 

+ f{Aq^q, + X, Aq,q, + X, 2v{x) + v{-X - Aq^q,), 2v{x) + v{-X - Aq^q,))] , 



(C24) 
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X, Ag^g^ - X, V{x) + 2v{Aq^q, - x),2y{x) + v{Aq^q, ~ x)) 
X, Aq^q^ - X, 2iy{x) + v{Aq^q> - x),2u{x) + v{Aq^q, - x))] 
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